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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 

TECHNICAL  NOTE  2797  ' 

A  STUDY  OF  THE  TRANSIENT  BEHAVIOR  OF  SHOCK  WAVES 
IN  TRANSONIC  CHANNEL  FLOWS 
By  Robert  V.  Hess 

SUMMARY  ♦  ' 


The  accuracy  of  the  result  obtained  in  a  fundamental  paper  by 
Kantrovitz  (NACA  TN  1225)  that  a  small  short-time  lowering  of  the  back 
pressure  in  steady,  shock-free,  transonic  diffuser  flow  causes  a  sta¬ 
tionary  or  trapped  shock  to  form  near  the  critical  sonic  channel  throat 
is  investigated  by  considering  the  contribution  of  a  higher-order  term 
in  the  short-time  calculations  which  was  neglected  in  Kantrowitz*s 
paper.  In  this  more  accurate  approximation  to  the  short-time  effects, 
the  shock  is  no  longer  stationary  or  trapped  unless  it  is  supported  by 
a  negative  steady-flow  back  pressure.  The  inclusion  of  the  higher-order 
term  in  the  short-time  calculations  avoids  the  use  of  approximate  quasi- 
steady-flow  considerations  for  the  complete  diffuser  flow  to  increase 
the  accuracy  of  the  shock  motion,  as  was  required  in  Kantrowitz’s  paper. 
In  a  broad  sense,  the  present  paper  offers  a  firmer  basis  for  the  short- 
time  approach  originated  in  Kantrowitz*s  paper. 

The  present  results  transform  into  those  previously  reported  in  NACA 
TN  1878  for  amplitudes  that  are  small  con^ared  to  the  difference  between 
local  and  critical  sonic  velocities  of  the  channel  flow. 


INTRODUCTION 


In  Kantrowitz^s  paper  (ref.  l),  the  time-dependent  shock  behavior 
produced  by  lowering  the  back  pressure' and  the  stability  of  steady 
shock-free  transonic  diffuser  flows  are  treated  by  dividing  the  time 
history  of  the  phenomena  into  short-time  effects  and  long-time  effects. 
The  short-time  effects  are  analyzed  for  a  single  distinct  upstream  dis¬ 
turbance  or  short  pulse  for  which  the  equations  describing  the  character 
of  unsteady  flow  are  greatly  simplified.  For  additional  simplification 
of  the  short -pulse  equations  the  highest -order  term  is  neglected.  The 
long-time  effects  are  concerned  with  transitory  phenomena  which  occur 
in  the  interval  between  the  end  of  short-time  phenomena  and  the  final 
steady  flow  state  in  the  channel.  The  calculations  indicate  that  the 
short-time  effect  of  a  small-amplitude  expansion  pulse,  produced  by  a 
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small  short-time  lowering  of  the  hack  pressure,  consists  of  the  forma¬ 
tion  of  a  shock  from  the  crest  of  the  expansion  pulse  which  becomes  sta¬ 
tionary  or  trapped.  In  an  originally  shock-free  flow,  the  occurrence  of 
a  stationary  or  trapped  shock  unsupported  by  a  negative  back  pressure  is 
in  disagreement  with  steady- flow  solutions  for  stationary  shocks,  and  the 
accuracy  of  shock-velocity  considerations  has  therefore  to  be  increased. 
The  more  accurate  shock  calculations  are  made  in  reference  1  by  means  of 
an  approximate  application  of  the  involved  long-time  effects  and  the 
steady-flow  back  pressure  at  the  end  of  the  diffuser,  with  the  aid  of  a 
convenient  quasi-steady-flow  approach.  The  result,  known  from  experi¬ 
ments,  is  that  the  originally  stationary  or  trapped  shock  will  consume 
the* short-time  expansion  pulse  if  the  back  pressure  is  that  of  the  shock- 
free  steady  diffuser  flow  and  that  the  shock  has  to  move  to  a  new  posi¬ 
tion  if  the  back  pressure  is  reduced. 

In  the  present  paper  the  increase  in  accuracy  of  the  shock-velocity 
calculations  is  made  in  a  different  way.  Since  the  result  of  the  sta¬ 
tionary  shock  in  reference  1  is  obtained  by  neglecting  a  higher-order 
term  in  the  equations  for  short -time -pulse  motion,  the  contribution  of 
this  term  to  the  order  of  accuracy  of  the  shock  motion  should  be  inves¬ 
tigated  before  the  approximated  contributions  due  to  the  basically  dif¬ 
ferent  quasi-steady-flow  considerations  are  taken  into  account.  Such  an 
approach  shows  whether  the  short-time  considerations  directly  result  in 
a  stationary  or  trapped  shock  without  application  of  a  negative  back 
pressure.  * 


SYMBOLS 


a 

u 

X 

A 

^ulse 

t 


velocity  of  sound 
flow  velocity 

distance  along  channel  axis 
cross-sectional  area  of  channel 
pulse  area 

time 

parameters  of  characteristic  families)  used  as  quantities  for 
measuring  amplitudes  of  unsteady-flow  disturbances.  Also 
used  as  mere  labels  for  distinction  between  downstre^  and 
upstream  disturbances.  The  expressions  "disturbance"  and 
"pulse"  are  used  interchangeably  in  the  present  paper. 


P,  Q 
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a*  critical  sonic  channel  flov  velocity 

M  Mach  niomber,  u/a 


m  =  -  1 

Xs  location  of  initial  shock  formation 


^LE  coordinate  of  leading  edge  of  pulse 

Xipg  coordinate  of  trailing  edge  of  pulse 

Vi^  velocity  immediately  ahead  of  trailing  shock 

U2  velocity  immediately  "behind  trailing  shock 

7  ratio  of  specific  heats 

p  density 

*p  pressure 

Subscript: 

o  steady-flow  values 

A  prime  designates  deviations  from  steady  channel  flow. 


ANALYSIS 

Background  for  Present  Analysis 


As  a  time -dependent ,  upstream-moving  disturbance  or  pulse  moves 
through  the  steady  flow  gradient  in  a  diffuser  that  is  free  of  other 
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disturbances,  the  shape  of  the  pulse  is  changed  and  reflected  vaves  are 
produced  inside  the  pulse.  The  growth  in  amplitude  of  the  reflected 
downstream  waves  from  zero  magnitude  at  the  leading  edge  of  the  upstream 
pulse  to  the  value  when  leaving  the  pulse  depends  on  the  magnitude  of 
the  flow  gradient  and  the  time  the  reflected  waves  are  permitted  to  grow 
inside  the  pulse  or  on  the  length  of  the  pulse  (see  also  refs.  1  and  2). 

•  For- a  given  steady  flow  gradient  a  certain  pulse  length  can  always  be 
found  in. which  the  values  of  the  reflected  amplitudes  grown  from  zero 
inside  the  pulse  are  negligible  conpared  to  the  original  pulse  an5)litude. 
The  original  pulse  and  its  reflected  pulse  that  is  built  up  over  this 
length  and  leaves  the  original  pulse  can  thus  be  treated  separately. 

This  separability  of  original  and  reflected  pulses  holds  for  pulses 
which  are  locally  superposed  on  the  steady  flow  gradient  (the  extent 
of  the  local  neighborhood  depends  on  the  magnitude  of  the  steady  flow 
gradient ) . 

Disturbances  of  this  nature  have  the  great  advantage  of  presenting 
a  case  for  which  the  unwieldy  equations  for  unsteady  flow  disturbances  can 
be  solved  with  relative  simplicity.  The  application  of  the  localized 
pulses  is  mainly  in  the  study  of  disturbances  that  influence  the  steady 
channel  flow  for  only  a  comparatively  short  time.  In  contrast,  the  long¬ 
time  effects  which  lead  to  the  final  steady  flow  state  depend  on  the 
Influence  of  the  complete  diffuser  flow  on  the  reflections  and  repeated 
reflections  from  the  original  pulse.  Pulses  which  can  be  treated  inde¬ 
pendently  of  their  reflections  are  designated  as  short  in  the  subsequent 
development. 

The  basic  equations  for  motion  of  isentropic  unsteady  flow  dis- 
^  turbances  are  conveniently  expressed  in  terms  of  the  characteristic 

parameters  P  =  u  +  • — —  a  and  Q  =  u  -  y  a  (see  refs.  1  and  2)  as 

1-^  =  0  (1) 

dt  dx 


iS  _  ua  =  0  '  (2) 

dt  dx 

where 

dP  3P  /  ^ 

^  +  (u  +  a) — 

dt  St  Sx  ■  , 
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and 


The  velocities  of  propagation  of  P  and  Q  are  given,  respectively,  by 

dx 


dt 


=  u  +  a 


(3) 


and 


to 

dt 


u 


a 


(4) 


Since  the  unsteady  pulses  considered  represent  deviations  from  steady 
flov,  the  quantities  u  and  a  are  expressed  in  the  form 


¥ 


U  =  Uq  +  u' 


a  =  a^  +  a' 


and 


Q'  =  u'  - 


The  characteristic  parameters  P'  =  u*  +  ^ 

2  ,  ,  7-1 

a  ,  based  on  deviations  from  steady  flow,  which  are 


77! 

associated  with  downstream  and  upstream  motions,  respectively,  in  the 

channel  flow  are  conveniently  cbnsidered  as  the  ampli- 
downstream  and  upstream  pulses.  Equation  (1^)  for  the  potion  of 
^  ^  ^  pulses  indicates  that  the  foot  of  an  upstream  pulse  will 

travel  with  the  velocity  u^  -  a^,  whereas  a  given  amplitude  Q-  travels 
with  the  velocity  u  -  a  =  (uq  -  ao)  +  (u'  -  a‘). 

amplitude^^Q^^^are  motion  of  upstream  pulses  with  the 


dQ* 

dt 


r- 

"(1  Q')  + 


-  Q') 


dx 


(5) 
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+  1-1-2  P' 


+  2L1^Q- 

4 


These  equations  are  subsequently  simplified  ^y  the  ° 

short-pulse  concept;  that  is,  the  amplitudes  P  of 

inside  the  pulse  are  considered  to  he  neglxgible^compared  to  the  ai|5.1i- 
tude  Q’  of  the  pulse  under  consideration.  Setting  ^  ^  . 

is  equivalent  to  letting  u'  =  -5a';  therefore,  Q  -  ”+pvTnc!  of  the 
The  amplitude  of  the  short  pulse  can  thus  e 

velocity  deviation  u*.  For  the  sake  of  mathematical  simplicity,  equa 
tion  (5)  is  further  reduced  for  use  near  Mq  =  1  t»y  assuming  that  th 
anrolitude  O'  or  2u'  is  small  compared  with  the  critical  sonic  velo- 
city.  (This  assumption,  however,  is  not  a  requirement  for  short  pulses.) 
Equations  (5)  and  (6)  are  thus  reduced  to 


du^  - 
dt  [1 


ii  +  yk 


■3^ 


=  Ur^  -  a, 


7  +  1  , 

+  - -  u' 


In  spite  of  the  simplification  of  the  equations  for  pulse  motion 
achieved  through  the  short-pulse  concept,  the 
still  too  complicated  to  permit  simple  solutions.  It  w 
enSe  1  tLt,^n  the  neighborhood  of  =  1,  the  higher-order  tem  in 

equation  (7),  (Mo^  -  neglected  compared  to 

(l  +  Mo)u’  ^  2  since  -  1  is  negligibly  small.  For  the  solution 

of  equations  (7)  and  (8),  furthermore,  a 

quantities  Uo,  ao,  and  x  has  to  be  known.  In  the  neighborhood  of 
M,  =1,  uo  -  ao  is  conveniently  expressed  in  terms  of  the  distance  x 

from  the  critical  sonic  section  of  the  channel.  In  re  ’  terms 

Taylor  series  development  of  Uq  -  Sq  in  x  is  broken  off  with  terms 

of  the  order  x.  The  development  of  Uo(x)  broken  off- with  the  term  o 

the  order  x  is 


«  /  \ 

UqCx)  =  Uo(0)  +  xf-^J 


„  a  *  +  — 

"  hi 
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vhere  =  & 


is  a  negative  quantity  following  the  notation  of 


reference  1,  where  x  is  positive  in  the  direction  of  decreasing  steady 
flow  velocity  Uq  (see  also  fig.  l). 

The  quantity  ao  is  related  to  Uq  hy  the  Bernoulli  equation 

Uo^  +  y4i:  (10) 

Development  of  the  square-root  expression  for  a.Q  and  subtraction 
from  Uq  yields 


Uo  -  ao  = 


7+1  X 
2  ^ 


Substitution  of  equation  (ll)  into  equations  (7)  and  (8)  without  the 
term  -  i)u'^^^ - 7  -  ij  results  in 


7+1  u 
2  b 


^  =  ZJlI  JL  +  21±J:  u' 

dt  2  hi  2 


The  solution  of  this  system  of  equations  is 


u'  =  Constant  (l4) 

The  pulse  distortion  near  Mq  =  1  thus  follows  a  family  of  hyperbolas, 

in  which  one  asymptote  has  the  slope  — -  of  the  velocity  gradient  and 

dx 

the  other  is  symmetric  to  it  with  respect  to  the  x-axis.  Figure  1  shows 
the  distortion  of  an  expansion  pulse  in  the  velocity  plane;  an  expansion 
pulse  is  chosen  since  its  shock  formation  is  of^^^pecial  interest.  A 
pulse  with,  for  example,  the  original  shape  ABCD  is  distorted  as  it 
approaches  Mq  =  1  and  takes  the  shape  ABCD.  In  the  course  of  this 

distortion  along  the  family  of  hyperbolas,  the  part  of  the  pulse  to  the 
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right  of  the  line  Xg  =  Constant  will  come  to  lie  at  a  larger  distance  x 
from  the  minimum  section  of  the  chainnel  than  other  parts  of  the  pulse. 

The  overhanging  part  CEF  of  the  pulse  cannot  physically  exist  since, 

when  the  crest  of  the  pulse  reaches  point  its  slope  is  infinite 

and  shock  fonnation  therefore  occurs.  A  triangular  pulse  with  a  trailing 
shock  is  finally  formed^  the  ‘leading  expansion  phase  of  the  pulse  is  at 
dUo 

the  slope  -  - - . 

du  * 

Because  of  the  occurrence  of  an  infinite  slope  desirable 

to  introduce  combinations  of  u’  and  x  which  will  not  undergo  such 
extreme  changes  as  Mq  =  1  is  approached.  .  The  pulse  area  is  such  a 
quantity.  As  is  indicated  in  reference  1,  the  pulse  area  is  sufficient 
to  determine  the  position  of  the  shock  trailing  a  triangular  pulse  near 
Mq  =  1,  because  the  shape  and  movement  of  the  leading  expansion  phase  as 

it  approaches'  Mq  =  1  are  known  and  simple.  For' the  case  discussed  in 
reference  1  ^that  is,  the  term  neglected^ 

the  pulse-area  growth  with  time  is  obtained  from  equations  ( 12 )  and  (13)- 
The  pulse-area  growth  with  time  is 


^(u'dx) 


where  the  differentiation 
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The  first  tvo  terms  cancel.  The  integral  of  the  remaining  term  from 
the  leading  edge  x^e  to  the  trailing  edge  of  the  pulse  is 

zero  since  the  corresponding  values  of  u*  are  equal  (zero)  at  these 
points.  The  pulse-area  growth  with  time  when  the  higher-order  term 


in  reference  1. 


is  neglected  is  thus  zero,  the  value  obtained 


Present  Analysis 


Mathematical  development.-  In  order  to  analyze  the  effect  of  the 

higher^order  term  (Mq^  -  pulse-area  growth 

near  the  critical  sonic  channel  section,  Uq  -  Eq  or  Mq  -  1  is 

developed  in  terms  of  the  distance  x  from  the  sonic  throat.  Since 
Mq  -  1  is  of  the  order  of  the  distance  x  made  nondimens ional  and  the 
distance  x  and  the  amplitude  u*  made  nondimens ional  are  ^assumed  to 
be  small  and  of  the  same  order,  the  lowest-order  contribution  of  the 


term  (Mq^  -  i^uY-^ - ^  is  of  the  order  x^.  The  condition  that 

\2Mq  4  / 

u*  and  X  be  of  equal  order  is  well-satisfied  for  a  triangular  pulse 
bounded  by  a  shock  near  Mq  ==  1  (u*  is  proportional  to  x  for  such  a 
pulse);  the  motion  of  such  a  triangular  shock-bounded  pulse  is  of  main 
interest  in  this  paper.  For  the  present  calculations  the  Taylor  series 

of  Uo(x),  broken  off  with  terms  of  the  order  x^  rather  than  x  as 
was  done  in  reference  1,  is  used: 


% 


Uo(x) 


(15) 


In  analogy  with  reference  1,  the  coefficients  are 


'  x=0  ^2 


and 


10 


MCA  TN  2797 


With  the  aid  of  the  Bernoulli  equation  (eq.  (lO)),  “  ^o 

expressed  in  terms  of  x  as 

* 


u 


o 


ao 


7+lx  .7  +  11  ^2 
"  2  5”  bg 


+ 


(7  -  1)2 

i 

8 

ao*bi2 

X  +  .  .  . 


(i6) 


In  order  to  obtain  Mq  -  1,  equation  (l6)  is  divided  by  aQ.  If,  for 
the  neighborhood  of  Mq  =  equation  (7)  is  written  in  the  modified  form 


u'  +  2(Mo  -  l) 


the  substitution  of  the  expressions  for  Uq  -  Uq  and  Mq  -  1  into  the 
modified  equation  (7)  and  equation  (8)  yields 


du’ 

7  +  1  u' 

7  +  1  u  *x 

(7  +  1)2 

+  il- 

-  7)(7  +  1) 

u’x 

dt 

CVl 

1^ 

OJ 

1 

8 

4 

bi^ao* 

and 

dx 

7  +  1  X 

.  7  +  1  X^  ,  ' 

7  -  1  ,  (7 

-  1)2' 

+I 

+  1  1 

u  * 

dt 

II 

ro 

4  b2 

8 

a^^bi^ 

2 

(17) 


(18) 


The  pulse -are a  growth 


n 


♦ 
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is  obtained  by  using  equations  (l7)  and  (18).  Wow 


^(u'dx)  =  dx  +  u'  —  dx 
dt  dt  dt 


where 


Go  Go 

^u'dx  -  ^xu'dx  + 
bi  bg 


ao*bi‘ 


xu'dx  + 


or 


C  C  2C 

^  u'dx  +  -2.  xu'dx  +  - xu’dx  +  Cou'du' 

^1  ^2  ao*bi2  3 


C  +  2C 

-^(u'dx)=  -i - ^  xu'dx  +  CoU'du' 

dt''  /  „  2  3 


ao*bi- 


(19) 


Cl  = 


(7  +  1)^  (3  -  r)(7  +  1) 

-  +  - 

8  4 


c„  .  2^  +  (?  - 

2  4  8 


and 


r  _  7  +  1 

Since  the  term  C^u^du’  disappears  vhen  the  integration  is  performed 
from  the  leading  edge  xj.-r  to  the  trailing  edge  of  the  pulse,  the 

pulse-area  growth  with  time  is 


+  2C2 


ao*h 


2 

1 


XU  *dx 


(20) 
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Now  the  pulse-area  growth  with  time  is  calculated  near  Mq  =  1, 
where  the  pulse  has  assumed  triangular  form.  The  reason  the  pulse  can 
still  be  considered  in  triangular  form  for  the  present  case  in  which 
terms  of  higher  order  than  those  in  reference  1  are  considered  for  the 
pulse-area  growth  is  as  follows:  In  the  present  approximation  the 
asymptotic  velocity  distribution  of  supersonic  steady  flow  approached 
by  the  leading  edge  of  the  expansion  pulse  is  no  longer  a  straight  line 
symmetric  to  the  linear  steady-flow  subsonic  velocity  distribution^ 
because  the  relation  between  velocity  and  cross-sectional  area  (or 
distance  x)  for  supersonic  and  subsonic  steady  isentropic  channel  flow 
is  asymmetric  with  respect  to  critical  sonic  velocity  (see  fig,  2  and 
the  appendix).  For  the  pulse  that  is  cut  off  by  the  discontinuous  shock 
at  a  small  value  of  x,  the  pulse  area  added  due  to  the  slight  deviation 
of  the  pulse  amplitude  u’  from  symmetric  conditions  will  be  of  higher 
order  than  that  being  considered  (see  fig,  2).  The  pulse-area  addition 
due  to  the  asymmetric  condition  of  the  present  approximation  can  thus  be 
neglected  and  the  pulse  can  be  considered  in  triangular  form.  The  integra¬ 
tion  in  equation  (20)  over  the  triangular  shape  is  simplified  by  the  fact 
that  the  expansion  phase  is  bounded  by  the  shock  which  coincides  with  the 
ordinate  u^^  ^ .  The  integration  therefore  only  has  to  be  performed  over 

the  leading  expansion  phase  from  x^  to  In  the  neighborhood  of 

Mq  =  1  the  leading  expansion  phase  can  be  considered  as  parallel  to  the 
asymptote  u'  =  -  ^  x  for  the  present  case  and  is  given  by 

u*  =  -  •  ^Le)-^  where  b^^  is  a  negative  quantity.  Substitution  of 

this  expression  for  u^  into  equation  (20)  yields  for  the  triangular- 
pulse  -area  growth  with  time 


d  ^TE 
—  u^dx  = 


2  +  2C2 

^1 


(x  - 


'  O  FlE 

L 


2(Ci  +  2C2)  (xjE  -  (xjE  ~  ^Le) 


As  Mq  - >  1  or  Xl2 


0  equation' (2l)  becomes 


d  te 


2(Ci+2C2)' 

u'dx  - - i - - —  x^ 

3ao*bi3 
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Since  b-|^  is  a  negative  quantity,  equation  (22)  shows  that,  as  the 
triangular  pulse  approaches  Mq  =  1,  its  area  grows  and  is  proportional 
to  the  cube  of  the  distance  from  the  critical  sonic  section.  Since  the 
area  of  a  triangular  pulse  Ap^j^gg  is 


A 

pulse 


2 


(23) 


equation  (22)  can  also  be  stated  in  the  form,  as 

■  1  '^(Vlse)  2  ^2 

^pulse  3 


(24) 


vhich  indicates  that  the  logarithmic  growth  rate  of  the  pulse  area  is 
proportional  to  the  distance  of  the  pulse-area  center  from  the  critical 
sonic  section. 


The  speed  of  the  trailing  shock  is  obtained  by  comparing  the  pulse- 
area  growth  from  equation  (21)  with  the  growth  obtained  by  differentiating, 
equation  (23)  with  respect  to  time 


(^ulse)  2  +  2C2 

(xte  -  Xle)^  ^  (x|rE  - 

^1  ao*bi2 

±jrj  2  3 

dxT-p  dx__ 

where  the  term  — ; —  still  has  to  be  determined.  The  speed  — —  of 
dt  dt 

the  foot  u'  =  0  of  the  leading  phase  of  the  expansion  pulse  (see 
,  eq.  (18))  is 


^LE 


7+11 
2  bi 


Xjp  +  Constant  x^^ 


2 


NACA  TN  2797 


14 


As  Mq - >1  or  Xyp, 

the  pulse  therefore  becomes 


dx 


- ^rp-g 

■>0,  the  speed  of  the  trailing  edge  of 


dx^g 

dt 


Cl  +  2C2 


Xte‘ 


(26) 


2 

In  terms  of  the  pulse  amplitude  u^'  =  -  —  the  speed  of  the 

trailing  edge  of  the  pulse  is 


dx™  Cl  +  2C2  ui 


.2 


dt 


12 


(7  +  1)(3  -  7) 
96 


,2 

IT 


(27) 


According  to  equation  (27),  the  trailing  edge  moves  away  from  the 
critical  sonic  channel  section  in  agreement  with  the  pulse-area  growth 
given  by  equations  (22)  and  {2h) , 

Correction  to  increase  accuracy  of  shock  velocity. -  In  the  develop¬ 
ment  of  equation  (27),  the  motion  of  the  discontinuous  trailing  edge  has 
been  identified  with  that  of  the  trailing  shocks  without  further  dis¬ 
cussion  of  the  process  in  which  the  shock  is  obtained  from  the  over¬ 
hanging  pulse.  In  making  no  special  issue  concerning  the  effect  of 
shock  formation  on  the  pulse  growth^  it  is  tacitly  taken  for  granted 
that  it  does  not  affect  the  magnitude  of  the  pulse  area  as  a  function 
of  time.  In  that  case  the  shock  is  obtained  by  averaging  the  pulse  area 
such  that  the  overhanging  area  CGM  is  equal  to  the  new  area  MEN  added 
by  the  shock  (see  fig.  l) .  This  averaging  process  for  obtaining  shocks 
is  known  to  apply  to  disturbances  of  small  amplitudes  u^  moving  through 
constant  cross  section  (ref.  3)«  The  applicability  of  this  averaging 
process  to  the  present  case  of  varying  cross  section  is  directly  con¬ 
nected  with  the  fact  that  the  pulse  will  remain  a  shock-bounded  triangle 
as  it  approaches  Mq  =  1.  More  specifically,  any  time  the  pulse  begins 
to  overhang  it  is  at  once  averaged  by  a  shock  (fig.  1  represents  an 
exaggerated  picture  of  the  process  of  shock  formation)^  that  is,  the 
averaging  process  by  the  shock  occurs  in  neighboring  sections.  The 
shock-wave  influence  can  thus  be  treated  as  independent .of  variation  of 
cross  section  (the  exact  proof  lies  in  the  fact  that  the  thickness  of 
the  shock  is  zero  and  its  pressures  are  bounded) . 

Since  an  averaging  process  applies  to  the  small  amplitudes  u’,  and 

2 

equation  (27)  for  the  shock  speed  contains  a  term  of  the  order  u^^*  , 
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the  problem  arises  whether  the  averaging  process  which  gives  the  .first- 
order  effect  of  the  shock  has  to  be  corrected  with  a  higher-order  term. 
It  is  indicated  in  the  preceding  paragraph  that  the  process  of  shock 
formation  can  be  treated  independently  of  the  variation  in  cross  section 
thus,  the  higher-order  effects  in  the  expressions  for  shock  speed  in 
constant  cross  section  can  be  used.  The  development  for  a  downstream 
shock  is  (in  the  notation  of  ref.  3,  eq.  (72. 05)) 

U  =  Uo  +  Co  +  -  Uo)  +  (28) 

Also,  it  can  be  found  from  reference  3,  equations  (72.06) 
and  (72.03),  that  the  sum  of  the  terms  Uq  +  Cq  +  ^(u  -  Uq)  is 

equivalent  to  the  average  velocity  ^^Uq  +  +  u  +  c).  In  other  words, 

the  first  three  terms  in  equation  (28)  deal  with  the  averaging  of  the 
speed  of  the  overhanging  pulse  by  that  of  the  shock.  ^Note  that,  since 
in  the  present  case  upstream  waves  are  being  considered,  the  averaged 

shock  velocity  is  i(uo  -  Cq  +  u  -  c).) 


Introduction  of  the  higher-order  effect 


(7  +  1)^  ^1'^ 
32  ao* 


(in  the 


present  notation)  effectively  increases  the  strength  and  speed  of  the 
trailing  shock  which  moves  upstream  during  its  formation  and  thus  has 
to  be  subtracted  from  the  shock  speed  based  on  averaging  the  pulse-area 
growth  (eq.  (27)).  The  proper  speed  of  the  trailing  shock  thus  is 


dt 


y  +  i[3  -  7 
32  L  3 


-  (7  +  1) 


-IJ  — 
J  ao* 


7+1  ^l'^ 

2k  ao* 


(29) 


The  negative  speed  indicates  that  the  shock  tends  to  consume  the  short- 
time  expansion  pulse  as  the  pulse  approaches  =  1.  Note  that  the 
ratio  of  the  proper  negative  shock  speed  (eq.  (29))  and  the  averaged 

shock  speed  (eq.  (27))  is  ^  =  -3.5  (for  air,  7  =  1.4).  The  proper 


l6 
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negative  shock  speed  is  thus  several  times  larger  than  the  averaged 
positive  speed  based  on  the  pulse-area  grovth. 

Conditions  for  stationary  trailing  shock. -  An  independent  check  of 
sign  and  magnitude  of  the  present  results  is  simple  to  make  with  the 
aid  of  steady- flow  considerations  for  a  stationary  trailing  shock.  In 
order  to  bring  the  pulse -consuming  shock  to  rest,  an  upstream- moving 
pulse  area  has  to  be  fed  into  the  shrinking  pulse  area  bounded  by  the 
moving  shock.  The  amplitude  of  this  fed-in  pulse  area  (see  fig.  2) 

has  to  be  equal  to  the  back  velocity  directly  behind  the  stationary  shock 
obtained  from  steady-flow  considerations  in  the  region  extending  from 
the  critical  sonic  throat  to  the  rear  of  the  shock  (see  the  appendix). 

The  stationary  shock  used  for  such  a  check  does  not  represent  the  actual 
steady-flow  solution  that  has  to  be  adjusted  to  a  certain  back  pressure 
at  the  end  of  the  complete  diffuser;  however,  it  represents  a  possible 
steady-flow  solution  for  a  different  back  pressure  at  the  end  of  the 
complete  diffuser.  The  magnitude  of  the  error  in  the  actual  steady-flow 
back  pressure  (velocity)  at  the  end  of  the  diffuser  depends  on  the 
importance  of  the  influence  of  the  long-time  effects ,  governed  by  the 
repeated  reflections,  leading  to  the  steady  flow  state  in  the' complete 
diffuser.  The  present  check  of  the  short-time  results  gains  in  physical 
significance  because,  as  discussed  in  a  subsequent  section,  the  long¬ 
time  effects  will  exert  relatively  minor  influence. 

The  back  velocity  necessary  to  bring  the  shock  motion  to  rest  is 
found  from  the  equivalent  condition  that  the  sum  of  the  pulse -area 
shrinkage  per  unit  time  due  to  the  shock  and  the  pulse-area  gain  per 
unit  time  due  to  application  of  the  (still  unknown)  positive  back  velo¬ 
city  immediately  behind  the  shock  is  zero.  The  amplitude  of  the 
reflected  downstream  pulse  is  shown  subsequently  to  be  negligible  com¬ 
pared  to  the  back  velocity.  The  pulse-area  shrinkage  per  unit  time  of 
the  triangular  pulse  is  given  for  Xj^ — ->0  (eq.  (25))  by 


'pulse 

dt 


™  dt 


2 

Introduction  of  equation  (29)  with  u^^*  ”  “  “  yields  for 

pulse-area  shrinkage  ^b2_  is  a  negative  quantity^ 


‘^(■^pulse)  1  (7  +  1)7  „  3 

1 i 
3  ■u_3fl  * 


C30) 
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The  rate  at  which  the  pulse  area  due  to  hack  velocity  enters  the  tri 
angular  pulse  at  the  shock  location  is  given  hy 


^  du' 


because  the  speed  at  which  the  various  parts  of  the  hack -velocity  pulse 
enter  the  shock  is  their  velocity  in  the  negative  x-direction.  Equa¬ 
tion  (18)  for  ^  is  then  applied,  and 
dt 


2  ^ 


27+  IpTE  \  ^  P 

— o — +  uV  +  Constant  Xm^'^ 


2  Vhi 


^2'  -  ^2'^  -  Constant  x^g% 


If  the  sum  of  the  triangular -pulse -area  shrinkage  (eq.  (30))  and  the 
pulse-area  addition  based  on  the  hack  velocity  Ug ’  is  set  equal  to 

zero,  where  x^  =  _  ^ 


,2(17  +  1)7 


2ha^* 


+  Constant  u. 


7+1 

4 


In  equation  (31)  the  back  velocity  up '  is  to  he  considered  negligibly 
small  compared  with  u^^'.  The  condition  of  smallness  of  Up  ’  was 

actually  introduced  previously  by  the  fact  that  the  shock  was  assumed 
to  cut  off  the  full  height  of  the  pulse  triangle.  (In  order  to  include 
the  possibility  of  values  of  Up’  which  are  of  the  order  of  u^',  the 

whole  process  of  shock  formation  would  have  to  be  considered,  during 
which  the  shock  would  cut  off  only  a  finite  part  of  the  pulse.  For  the 
present  problem  the  final  stages  of  pulse  shrinkage,  represented  by  the 
pulse-area  triangle,  are,  however,  of  primary  interest.)  Neglect  of  up 

compared  with  u^'  in  equation  (31)  yields  for  the  back  velocity  up ' 


6  a^* 
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It  was  previously  stated  that  the  amplitude  of  the  positive  hack 
velocity  U2'  that  will  bring  the  shock  to  rest  would  have  to  he  equal 
to  the  hack  velocity  required  immediately  behind  a  stationary  steady- 
flow  shock.  This  back  velocity  is  obtained  by  subtracting  the  velocity 
increment  through  a  steady  shock  near  Mq  =  1  from  the  increment  between 
supersonic  and  subsonic  velocities  corresponding  to  the  same  channel 
area  near  Mq  =  1  for  isentropic  flow.  In  the  appendix  this  difference 
is  shown  to  be  equal  to  the  value  in  equation  (32). 

This  calculation  of  the  back  velocity  U2'  is  especially  simple 
because  the  amplitude  of  the  reflected  downstream  pulse,  which  is  built^ 
up  in  the  triangular  upstream  pulse,  can  be  neglected.  Under  such  condi¬ 
tions  the  back  velocity  required  to  bring  the  shock  (trailing  the  tri¬ 
angular  pulse)  to  rest  does  not  have  to  be  superposed  on  the  reflected 
amplitude  but  can  be  directly  superposed  on  the  foot  of  the  trailing- 
shock  amplitude.  The  proof  for  the  permissibility  of  the  neglect  is  now 
given. 

The  arqslitude  AP'  of  the  built-up  reflected  short  pulse  near  Mach 
number  1  is  obtained  from  the  equation  for  a  downstream  or  P'  pulse. 

A  numerical  estimate  of  the  reflected  amplitude  is  given  in  reference  Ij 
however,  its  order  of  magnitude  is  not  specified  there.  The  derivation 
is  thus  repeated  here  with  different  emphasis.  The  growth  dP  of  a 
reflected  P'  wave  is  given  by  (for  air,  7  =  1-^) 


dP' 

dt 


5 


vhere  ^*20"  neglected.  For  consideration  of  as  the  growth 

of  the  reflection  of  a  qhort  Q'  pulse,  P'  is  neglected  compared  with 
Q';  thus,  near  Mq  =  1, 


dP' 

“dt 


For  a  downstream  pulse the  relation  between  t  and  x  is  given  by  the 
downstream  speed  of  the  pulse  foot  (for  the  present  order  of  accuracy) 


dx 

dt 


+  ^o 


as 
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Introducing  1  -  Mq  =  Constant  x  (eq,.  (l6))  to  obtain  the  first-order 
effect  for  the  reflected  amplitude  and  noting  that  Q'  =  2u'  yields 

du- 

dP’  =  Constant  u'x  - 2 — 

Uq  +  ao 


Near  =  1^  u^  +  a^  2aQ*  and  writing 


dUo  T 

dUo  =  dx  =  ^  dx 


results  in 


dP^  =  Constant  u’x  dx 

The  reflected  amplitude  is  obtained  by  usir^  the  relation 

u’  =  Constant ^x  -  from  the  triangular  pulse ^  integrating  from 

Xle  to  Xfj^,  and  letting  Xj^ - >0 


AP’  =  Constant  U]_’^ 


The  reflected  amplitude  produced  by  the  entropy  change  in  the  shock 
is  also  proportional  to  the  third  power  of  the  velocity  change  in  the 
shock  (ref,  3) •  The  amplitude  of  the  reflected  pulse  is  thus  negligible 
compared  with  the  back  velocity  Up’  which  is  of  second  power  in  u-i  ' 


dUo 

’Tx 


2. 

b 


has  a  reasonable  value).  In 


(provided  the  flow  gradient 

\  U  A 

contrast,  as  previously  indicated,  the  pulse-area  growth  and  the  related 
ref lected-pulse  area  cannot  be  neglected  under  considerations  of  the  same 
order. 


It  should  be  pointed  out  that  the  equations  used  in  this  analysis 
and  in  reference  1  for  the  determination  of  the  downstream  reflection 
leaving  a  short  upstream  pulse  are  based  on  isentropic  considerations. 

They  therefore  do  not  apply  directly  to  the  pulse-area  shrinkage  due  to 
the  motion  of  the  proper  shock  but  rather  to  the  pulse -area  growth  averaged 
by  the  shock.  The  positive  speed  of  the  averaging  shock  (eq.  (27)) 
and  the  negative  speed  of  the  proper  shock  (eq.  (29))  are,  however/  of 
the  same  order  of  magnitude  (their  ratio  of  -3.5  is  not  sufficiently  large 
to  affect  the  order).  The  reflected  amplitude  AP’  based  on  motion  of 
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the  proper  shock  is  thus  of  the  same  order  of  magnitude  as  the  ampli- 
tude  AP'  based  on  pulse-area  averaging. 


DISCUSSION  OF  RESULTS 


The  result  has  been  obtained  that  the  consumption  of  the  expansion 
pulse  by  the  shock  does  not  have  to  wait  for  long-time  effects  depending 
on  repeated  reflections  and  leading  to  the  final  steady  flow  state  in 
the  complete  diffuser  but  that  the  consumption  begins  as  soon  as  the 
shock  has  formed.  The  speed  of  pulse  consumption  was  considered  for  the 
case  in  which  the  shock  is  already  fully  formed  (for  the  nature  of  the 
shock  formation,  see  fig.  l),  and  the  pulse  was  conveniently  considered 
to  be  in  the  immediate  neighborhood  of  Mach  number  1  where  the  leading 
expansion  phase  of  the  pulse  approaches  zero  speed.  The  fact,  however, 
that  the  proper  negative  speed  of  the  pulse-consuming  shock  (eq.  (29)) 
is  larger  by  the  appreciable  factor  of  3-5  than  the  shock  speed  in  the 
opposite  direction  based  on  averaging  the  positively  growii^  pulse  area 
indicates  that  the  pulse  consumption  begins  when  the  shock  is  still  in  the 
last  stages  of  formation  and  the  leading  expansion  phase  of  the  pulse 
has  a  small  but  finite  speed. 

The  long-time  effects  exert  their  influence  through  the  repeated 
reflections  from  the  complete  diffuser  flow,  which  in  turn  affect  the 
shock  velocity  directly  through  a  change  in  back  velocity  immediately  ^ 
behind  the  shock.  The  amplitudes  of  the  first  downstream  reflections 
have  been  shown  previously  to  be  small  compared  to  the  back  velocity 
near  Mach  number  1.  The  repeated  reflections,  which  are  the  basis  of 
“the  long-time  effects,  are  small  compared  to  the  first  downstream  reflec¬ 
tions  and  they  move  upstream  at  a  low  speed.  The  negligible  contribution 
of  the  long-time  effects  indicates  that  the  back  velocity  immediately 
behind  the  shock  obtained  from  short -time -pulse  considerations  is  essen¬ 
tially  equal  to  that  corresponding  to  the  actual  steady-flow  back  pres¬ 
sure  at  the  end  of  the  diffuser.  The  present  higher-order  approximation 
to  the  short -pulse  considerations  for  small  shocks  hear  Mach  number  1 
thus  Indicates  that  shock  motion  occurs  immediately  after  shock  formation 
and  is  essentially  the  same  as  that  due  to  the  actual  back  pressure  at 
the  end  of  the  diffuser.  The  approximate  considerations  of  the  long-time 
effects  and  of  the  actual  steady-flow  back  pressure  which  had  to  be  made 
in  reference  1  are  thus  avoided  in  the  present  paper. 

In  reference  1  a  procedure  that  is  somewhat  similar  to  that  in  the 
present  paper  is  used  in  that  the  channel  end  conditions  are  also  applied 
directly  behind  the  shock.  However,  in  view  of  the  fact  that  in  refer¬ 
ence  1  the  increased  accuracy  of  the  shock -velocity  calculations  is  not 
balanced  by  an  increased  accuracy  in  the  short -time -pulse  equations  (7) 
and  (8),  a  true  measure  of  the  negligible  smallness  of  the  long-time 
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effects  of  repeated  reflections  is  lost.  This  situation  is  the  actual 
reason  why  the  increased  accuracy  of  the  shock  velocity  in  reference  1 
.  has  to  be  gained  through  discussion  of  the  various  steady-flow  conditions 

at  the  end  of  the  diffuser  and  why  the  approximate  treatment  of  the  long¬ 
time  effects  has  to  be  considered  as  a  "simplifying  assumption"  and  the 
long-time  effects  still  cannot  be  truly  neglected.  It  should  be  empha¬ 
sized  that  in  the  present  short-time  considerations  the  speed  of  the 

shock  which  consumes  the  pulse  is  of  the  order  x^;  whereas  the  speed 
of  the  pulse  before  the  leading  expansion  phase  approaches  zero  speed 
is  of  the  order  of  x  or  u'  (eg.  (I8)).  The  present  results  thus 
agree  with  those  obtained  from  approximate  quasi-steady-flow  considera¬ 
tions  in  reference  1  in  that  the  speed  of  pulse  consumption  for  small 
shocks  will  be  small  compared  with  the  speed  of  pulse  approach.  The 
present  paper,  however,  indicates  that  a  stationary  shock  does  not  occur 
unless  a  negative  back  pressure  is  applied  and  that  the  shock  begins  to 
consiame  the  pulse  immediately  after  the  shock  is  formed.  This  paper  , 
also  avoids  the  introduction  of  the  approximate  quasi-steady  considera¬ 
tions  (dealing  with  the  complete  diffuser  flow)  for  proof  of  the  condi¬ 
tions  of  pulse  consumption  and  thus  gives  a  firmer  basis  for  the  funda¬ 
mental  short-time  approach  in  reference  1.  The  fact  that  for  the  diffuser 
flow  the  shock  is  not  trapped  after  formation  but  tends  to  consume  the 
expansion  pulse  directly  after  formation  has  bearing  on  the  problem  of 
existence  of  shocks  in  transonic  flows  about  bodies. 

•  .  The  conclusions  concerning  the  nature  of  shock  formation  do  not 

affect  the  important  result  obtained  in  reference  1  that  short  compres¬ 
sion  pulses  moving  upstream  can  be  compensated  by  the  area  of  a  sta- 
’  tionary  but  still  short  (as  defined  herein)  expansion  pulse  which  is 

part  of  the  undisturbed  steady  supersonic  flow  bounded  by  a  shock.  (For 
a  given  shock  amplitude,  an  increase  in  throat  length  increases  the  area 
of  the  compensating  expansion  pulse.) 

The  preceding  results  of  pulse-area  growth  are  now  compared  with 
those  in  reference  2.  For  purposes  of  mathematical  facility,  equa¬ 
tions  (7)  and  (8)  were  simplified  by  the  restriction  u'  «  Uq  -  a^. 

The  significance  of  this  simplification  is  quickly  seen  by  substituting 
equation  (8)  Into  equation  (7),  which  results  in  a  differential  equation 
for  u'  in  terms  of  Mq,  for  which  separation  of  variables  is  possible. 
The  restriction,  however,  has  the  following  effect  on  the  pulse  distor¬ 
tion  for  which  the  term  can  be  neglected) 

namely,  equations  (12)  and  (I3)  for  the  pulse  distortion  are  reduced  to 

du'  _  7  +  1  u' 

'  dt  “  ■  2  bi 

% 

dx  _  7  +  1  X 
dt  ■  2  b^ 


(33a) 
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Their  solution  is 

xu^  =  Constant  (33"b) 

For  such  a  distortion  the  pulse  will  not  overhang  (see  preceding  analysis 
as  it  approaches  Mq  =  1>  and  no  shock  formation  will  occur.  Note  that 
the  restriction  u’  «  x  causes  equations  (33)  ‘the  pulse  distortion 
to  break  down  as  x  — >0,  because  u'  — the  restriction  requires 
that  the  distance  of  the  pulse  from  the  sonic  channel  throat  remain  of 
the  order  x.  Athough  such  a  simplified  approach  will  give  the  proper 
order  of  magnitude  of  the  pulse-area  growth  near  Mq  ==  1  without  the 
effects  of  the  shocks  it  cannot  treat  the  final  shrinkage  of  the  pulse 
bounded  by  a  shock.  In  order  to  derive  an  expression  analogous  to  that 
appearing  in  reference  2  for  the  rate  of  logarithmic  growth  of  the  pulse 
area,  equation  (20)  is  divided  by  the  integrated  pulse  area.  As  the 
leading  and  trailing  edges  of  the  pulse  are  permitted  to  approach  each 
other,  the  resulting  quotient  of  integrals  is  the  location  of  the  center 
of  the  pulse  (or  x)  times  a  constant 


J^LE 


In  reference  2,  the  pulse-area  growth  with  time  near  Mq  =  1  is 
not  given^  however,  the  pulse-area  growth  with  Mach  number  is  determined. 
Equation  (3h)  is  reduced  to  the  result  in  reference  2  as  follows:  A 
relation  between  t  and  x  is  obtained  from  . 
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ta  _  7+1  1 
dt  2 


(35) 


Substitution  of  equation  (35)  into  equation  (3^)  gives  for  the  pulse-area 
growth  with  distance  ' 


d  log  I  u'dx 

>rLE _ 

dx 


Cl  +  2C- 


7+1 


ao*bi 


(36) 


Since 


equation  (36)  becomes 


-  _  ^1  ^^2 

log  u'dx  =  - : 

J^LE  --4—  ao 


d 


(37) 


Furthermore,  near  Mq  =  1, 


2 

7  +  1 


dMo 


and 


Cl  =  - 


(7  +  1)^  (3  -  7)(7  +  1) 

8  h 


Co  =  7-1  +  (7  -  1)^ 

k  8 


Equation  (37)  becomes 


d  log 


u'dx 


=  3-7 
2(7  +  1) 


(38) 
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In  reference  2^  the  expression  d  log(u'dx)/Mo  has  the  same  signif¬ 
icance  as  the  expression  d  log(Ju'dxydMo  of  equation  (38).  In  agree¬ 
ment  with  equation  (74)  of  reference  2,  the  pulse-6irea  growth  with  Mach 
number  near  M©  =  1  is  thus  equal  to  I/3  for  7  =1.4  (air).  The 

pulse-area  growth  is  zero  for  a  hypothetical  gas  with  7=3  (one  degree 
of  freedom) . 

A  simple  physical  interpretation  of  the  pulse-area  growth  can  be 
given  by  noting  that  for  small  disturbances  the  pulse  area  Ju'dx,  if 
multiplied  by  the  density  Pq,  has  the  dimension  of  a  momentum!  Since 
the  momentum  has  the  dimension  of  the  product  of  force  and  time>  the 
growth  rate  of  the  pulse  area  is  proportional  to  a  force.  In  the  pre¬ 
sent  case  of  channel  flow  the  force  is  represented  by  the  axial  component 
of  the  normal  force  on  the  channel  wall.  The  normal  force  is  given  by 
the  integral  of  the  excess  pressure  p'  over  the  length  of  the  pulse 
and  is  of  finite  value  (note  the  connection  of  the  finite-normal -force 
integral  with  the  pulse  area).  As  the  pulse  approaches  Mq  =  1^  the 
slope  of  the  channel  wall  becomes  zero  and  thus  the  axial  component  of 
the  finite  normal  force  becomes  zero.  It  is  indicated  that,  in  order 
that  the  pulse-area  growth  with  time  be  zero,  the  entire  pulse  would 
have  to  move  to  the  critical  sonic  channel  section.  Since^  of  course, 
for  the  present  triangular  pulse  only  the  leading  edge,  not  the  center 
of  the  pulse,  approaches  asymptotically  the  sonic  section,  the  rate  of 
pulse-area  growth  (exclusive  of  the  shock  motion)  remains  slightly 
positive. 

The  simple  physical  picture  using  momentum  considerations  can  be 
also  used  for  an  explanation  of  the  pulse-area  growth  with  distance. 

The  pulse-area  growth  with  distance  is  of  the  dimension  Force  X  Time 

Force  x  Time  Distance 

7 - r - — — For  a  pulse  (or  pulse  part)  that  reaches  Mo  =  1, 

(u  -  a)  X  Time  •  o  ^ 

the  force  is  zero,  and  the  pulse  velocity  u  -  a  also  becomes  zero. 
According  to  these  considerations,  the  pulse-area  growth  with  distance 
near  Mq  =  1  should  thus  be  of  the  order  of  unity. 


CONCLUSIONS 


A  study  of  the  effects  of  a  small  short-time  lowering  of  the  back 
pressure  in  steady  shock-free  transonic  diffuser  flow  by  means  of  a 
higher  approximation  than  in  NACA  TN  1225  yields  the  following  conclusions 
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1.  The  more  accurate  approximation  of  the  present  paper  to  the 
short-time  effects  shows  that  the  shock  is  no  longer  stationary  or 
trapped  in  the  diffuser  unless  it  is  supported  by  a  negative  steady- 
flow  back  pressure)  the  result  thus  is  no  longer  in  disagreement  with 
steady-flow  solutions  for  stationary  shocks.  * 

2.  The  present  short-time  calculations  avoid  the  use  of  approxi¬ 
mate  quasi-steady-flow  considerations  for  the  complete  diffuser  flow 
to  increase  the  accuracy  of  the  shock  motion,  as  was  required  in 
Kantrowitz’s  paper  (NACA  TN  1225).  The  fundamental  considerations  in 
Kantrowitz’s  paper  are  thus  put  on  a  firmer  basis. 

3.  For  short-time  pulses  with  amplitudes  that  are  restricted  to 
values  that  are  small  even  compared  with  the  difference  between  local 
and  critical  sonic  velocities  of  the  channel  flow,  the  present  results 
transform  into  those  previously  reported  in  NACA  TN  I878. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  July  10,  1952. 
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APPENDIX 

THE  BACK  VELOCITY  DUE  TO  STEADY-FLOW 
#  SHOCK  LOSSES  NEAR  M  =  1 


The  back  velocity  is  obtained  by  subtracting  the  velocity  increment 
through  a  steady  shock  neai"  M  =  1  from  the  difference  between  super¬ 
sonic  and  subsonic  velocities  corresponding  to  the  same  channel  area 
near  M  =  1  for  isentropic  flow. 

The  velocity  Increment  through  the  shock  is  expressed  in  terms  of 
(ref.  3) 

Up  = 

(for  notation  used^  see  fig.  2)  as 


^^shock 


(Al) 


where  the  subscript  1  refers  to  the  supersonic  side  of  the  shock  and 
the  subscript  2  to  the  subsonic  side.  If 


(A2) 


is  substituted  into  equation  (Al),  the  result  is 


Au 


•shock 


m^ 


^1  +  m^ 


mi  -  i  m^^  + 


(A3) 


The  isentropic  channel-area  equation  (where  A*  is  the  critical 
sonic  channel  area) 


\A*y  m2|7  +  1\  2  / 


— ,  7+1 
7-1 
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is  altered  by  making  the  substitution 

- 

(7  +  1)  -  (7  -  1)M^ 

and  by  using  equation  (A2)  for  .  The  isentropic  channel-area  equa- 
tlon  then  becomes 

(Af  .  (1  *  „)-r3 

Development  of  equation  (a4)  results  in 
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